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>• 1 Measurements of if) and if)' production at the CERN SPS are compared to predictions 
based on a hadronic model of charmonium suppression. Detailed information is presented 
to facilitate comparison to other analyses. Sensitivity of these conclusions to the model 
parameters is discussed. 

1. INTRODUCTION 

Has the quark-gluon plasma been discovered at the SPS? The NA50 collaboration 
reported a dramatic decrease in if) production in Pb+Pb collisions at 158 GeV per nucleon 
0. Specifically, M. Gonin presented at this meeting a striking 'threshold effect' in the 
■?/>-to-continuum ratio as a function of a calculated quantity, the mean path length of 
the if) through the nuclear medium, L, see fig. 1. This apparent threshold has sparked 
considerable excitement as it may signal deconfinement in the heavy Pb+Pb system 0. 

In this talk we report on work in ref. || comparing Pb results to predictions using a 
hadronic model of charmonium suppression. We first demonstrate that the behavior in 
the NA50 plot is not a threshold effect but, rather, reflects the approach to the geometrical 
limit of L as the symmetric Pb+Pb collisions become increasingly central. When plotted 
as a function of the measured neutral transverse energy Et as in fig. 2, the data varies 
smoothly as in S+U measurements in fig. 4b below ||-|§. The difference between S+U 
and Pb+Pb data lies strictly in the relative magnitude. To assess this magnitude, we 
compare if) and if)' data to expectations based on the hadronic comover model |4]||. The 
curve in fig. 1 represents our calculation using parameters fixed earlier in Ref. M. Our 
result is essentially the same as the Pb+Pb prediction in 

Our primary intention in this work is to demonstrate that there is no evidence for a 
strong discontinuity between pA, S+U and Pb+Pb data. Secondly, we show that our 
model predictions agree with the new Pb+Pb data. The consistency of our predictions 
can be judged by comparing our old pA and S+U calculations with more recent NA38 
and NA51 data. However, we must be cautious about taking model details too seriously. 
Regarding the first point, we quote Maurice Goldhaber: ". . . absence of evidence is not 
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Figure 1. NA50 comparison [1] of Pb+Pb and S+U^ ip + X as functions of the average 
path length L, see eq. (3). B is the ip — > branching ratio. 



evidence of absence." Moreover, we do not attempt to show that our comover interpre- 
tation of the data is unambiguous - this is certainly impossible at present. Finally, the 
significance of our present results must be weighted by the fact that all pA and AB data 
are preliminary and at different beam energies. Whether the newer data is indeed 'im- 
proved' has yet to be established. One should be particularly cautious in taking literally 
the values of model parameters, such as the comover density. 



2. NUCLEONS AND COMOVERS 

The hadronic contribution to the suppression arises from scattering of the nascent ip 
with produced particles - the comovers - and nucleons [[|-f| . To determine the suppression 
from nucleon absorption, we calculate the probability that a cc pair produced at a point 
(6, z) in a nucleus survives scattering with nucleons to form a ip. The standard [pj[l0| 
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Figure 2. ^-to-continuum ratio as a function of Et from [1] compared to the comover 
prediction [3]. 



result is 

dz p A (b, z)a^ N } (1) 

where pa is the nuclear density, b the impact parameter and a^N the absorption cross 
section for -?/>-nucleon interactions. One can estimate Sa ~ exp{— a^^PoLA}, where La 
is the path length traversed by the cc pair. On average, La cannot exceed the nuclear 
radius Ra, since the pair can be produced anywhere within the nucleus. 

Suppression can also be caused by scattering with mesons that happen to travel along 
with the cc pair (see refs. in H). The density of such comovers scales roughly with Et- 
The corresponding survival probability is 

S co = exp{- J dm a co v re \} } (2) 

where n is the comover density and r is the time in the ip rest frame. We write S co ~ 
exp{— /3E T }, where j3 depends on the scattering frequency, the formation time of the 
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comovers and the transverse size of the central region, Rt, cf. eq. (8). 

To understand the saturation of the Pb data with L in fig. 1, we apply the schematic 



approximation of Ref. [10] for the moment to write 



t y TJ oc (S A S B S co ) ~ e-W oi e-^, (3) 



where the brackets imply an average over the collision geometry for fixed Et and ct(Et) = 
da/dE T . The path length L = (La + L B ) and transverse size Rt depend on the collision 
geometry. The path length grows with Et, but most stop at the geometric limit Ra + Rb 
in central collisions. Explicit calculations show that nucleon absorption begins to saturate 
for b < Ra, where Ra is the smaller of the two nuclei, see fig. 6 below. On the other 
hand, particle production and, consequently, Et continue to grow after L saturates due, 
e.g., to fluctuations in the number and hardness of NN collisions. Equation (2) falls 
exponentially in the 'super-central' regime, b < Ra, because (3 is essentially fixed. 

Important for understanding the NA50 results is the fact that saturation occurs at 
larger values of the average impact parameter for Pb+Pb system compared to S+U, see 
fig. 6 below. This result reflects the striking geometric difference between symmetric and 
asymmetric nuclear collisions. Specifically, in the S+U system there is large region, i?s < 
b < Rs + R\j, where varying b changes the volume of the interaction region dramatically. 
In contrast, Ra = Rb i n Pb+Pb collisions and i?pb ~ 6.6 fm is large. We show below that 
typical values of b corresponding to the Et region covered by NA50 data do not go far 
into the region i?pb < b < 2i?pb where [The importance of b being near zero for saturation 
to occur was overemphasised in SG's oral presentation of this talk.] 

3. J/ip SUPPRESSION 

In fig. 1, we compare the Pb data to calculations of the ?/>-to-continuum ratio that 
incorporate nucleon and comover scattering. The contribution due to nucleon absorption 
indeed levels off for small values of b, as expected from eq. (3). Comover scattering 
accounts for the remaining suppression. 

These results are predictions obtained using the computer code of Ref. |JJ] with param- 
eters determined in Ref. ||. However, to confront the present NA50 analysis [0J, we must 
perform the following updates: 

• Calculate the continuum dimuon yield in the new mass range 2.9 < M < 4.5 GeV. 

• Adjust the Et scale to the pseudorapidity acceptance of the NA50 calorimeter, 
1.1 < rj < 2.3. 

The agreement in fig. 1 depends strongly on these adjustments. 

We now review the details of our calculations, highlighting the adjustments as we go. 
For collisions at a fixed b, the ^-production cross section is 

a$ B (b) = a$ N J d 2 sdzdz' p A {s,z)p B {b - s,z') S (4) 

where S = SaSbS co is the product of the survival probabilities in the projectile A, target 
B and comover matter. The continuum cross section is 

°$vr ( & ) = o-X-Jd 2 sdzdz> PA (s, z)p B (b - s, z'). (5) 
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The magnitude of eqs. (4), (5) and their ratio is fixed by the elementary cross sections a^ N 
and cr^ u - . We calculate a^ N using the phenomenologically-successful color evaporation 



model [II]. The continuum in the mass range used by NA50, 2.9 < M < 4.5 GeV 



is described by the Drell-Yan process |fLz|| . To confront NA50 and NA38 data in the 
appropriate kinematic regime, we compute these cross sections at leading order using 
GRV LO parton distributions with a charm iT-factor K c = 2.7 and a color evaporation 
coefficient = 2.54% from [HJ and a Drell-Yan i^-factor Kdy — 2.4. Using these cross 



sections to construct the normalized ratios in fig. 2 corresponds to the first update. 
To obtain Et dependent cross sections from eqs. (4) and (5), we write 

a AB (E T ) = J d 2 bP(E T , b)a AB (b). (6) 

The probability P(Et, b) that a collision at impact parameter b produces transverse energy 
Et is related to the minimum-bias distribution by 

a min (ET) = Jd 2 bP(E T ,b). (7) 



We_parametrize P(E T} b) = C exp{-(E T - E T ) 2 /2A} } where E T {b) = eAf(b), A (6) = 
ueE T (b), C(b) = (27rA(6)) _1 and Af(b) is the number of participants (see, e.g., Ref. |§). 
We take e and uo to be phenomenological calorimeter-dependent constants. 

We compare the minimum bias distributions for total hadronic Et calculated using 
eq. (7) for e = 1.3 GeV and u = 2.0 to NA35 S+S and NA49 Pb+Pb data JT3l. The 



agreement in fig. 3a builds our confidence that eq. (7) applies to the heavy Pb+Pb system. 
Figure 3b shows the distribution of neutral transverse energy calculated using eqs. (5) 




Figure 3. Transverse energy distributions from eqs. (7, 8). The S-Pb comparison (a) 
employs the same parameters. 
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and (6) to simulate the NA50 dimuon trigger. We take e = 0.35 GeV, u = 3.2, and 
a u+L- ~ 35.2 pb as appropriate for the dimuon-mass range 2.9 < M < 4.5 GeV. The E? 
distribution for S+U — > + X from NA38 was described || using e = 0.64 GeV and 

u = 3.2 - the change in e corresponds roughly to the shift in particle production when the 
pseudorapidity coverage is changed from 1.7 < r) < 4.1 (NA38) to 1.1 < r\ < 2.3 (NA50). 
Taking e = 0.35 GeV is the second update listed earlier. Though not shown here, we 
comment that the calculations agree with the NA50 data at high Et, but differ somewhat 
at low Et due to the efficiency of the active target, which is lower for more peripheral 
collisions. 

We now apply eqs. (1,2,4) and (5) to charmonium suppression in Pb+Pb collisions. To 
determine nucleon absorption, we used pA data to fix a^N ~ 4.8 mb in Ref. ||. This 
choice is in accord with the latest NA38 and NA51 pA data, see fig. 4a. To specify comover 
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Figure 4. (a) pA cross sections [1] in the NA50 acceptance and (b) NA38 S+U ratios from 
the '91 [6] and '92 [1] runs. '92 data are scaled to the '91 continuum. The dashed line 
indicates the suppression from nucleon absorption alone. 



scattering ||, we assumed that the dominant contribution to ip dissociation comes from 
exothermic hadronic reactions such asp + ^— > D + D. We further took the comovers 
to evolve from a formation time r ~ 2 fm to a freezeout time t> ~ -Rr/^rei following 
Bjorken scaling, where v ie \ ~ 0.6 is roughly the average if) — p relative velocity. The 
survival probability, eq. (2), is then 

S co = exp{-cr co w re in r ln( J R T /t; re iro)} (8) 

where a co ~ 2<j^at/3, Rt ~ Ra and n is the initial density of sufficiently massive p, uj and 
r] mesons. To account for the variation of density with Et, we take n = uqEt / 'Et(0) ||. 
A value n Q = 0.8 fm -3 was chosen to fit the central S+U datum. Since we fix the density 
in central collisions, this simple ansatz for S co may be inaccurate for peripheral collisions. 
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We expect the comover contribution to the suppression to increase in Pb+Pb relative 
to S+U for central collisions because both the initial density and lifetime of the system 
can increase. To be conservative, we assumed that Pb and S beams achieve the same 
mean initial density. Even so, the lifetime of the system essentially doubles in Pb+Pb 
because Rt ~ Ra increases to 6.6 fm from 3.6 fm in S+U. The increase in the comover 
contribution evident in comparing figs. 2 and 4b is described by the seemingly innocuous 
logarithm in eq. (8), which increases by ~ 60% in the larger Pb system. 

In Ref . || , we pointed out that comovers were necessary to explain S+U data from the 
NA38 1991 run 0. Data just released [|IJ from their 1992 run support this conclusion. 
The '91 if) data were presented as a ratio to the dimuon continuum in the low mass range 
1.7 < M < 2.7 GeV, where charm decays are an important source of dileptons. On the 
other hand, the '92 if) data are ratios to the Drell-Yan cross section in the range 
1.5 < M < 5.0 GeV. That cross section is extracted from the continuum by fixing the 
-fT-factor in the high mass region ||14|| . To compare our result from Ref. || to these 



data, we scale the '92 data by an empirical factor. This factor is ~ 10% larger than our 
calculated factor a^y (92)/<J^ (91) ~ 0.4; these values agree within the NA38 systematic 
errors. [NA50 similarly scaled the '92 data to the high-mass continuum to produce fig. 1.] 
Because our fit is driven by the highest Ej- datum, we see from fig. 4b that a fit to the 
'92 data would not appreciably change our result. Note that a uniform decrease of the 
ratio by 10% would increase the comover contribution needed to explain S+U collisions. 

NA50 and NA38 have also measured the total ^-production cross section in Pb+Pb Jl]] 
and S+U reactions J7J. To compare to that data, we integrate eqs. (4, 6) to obtain the total 
(cr/AB)^ = 0.97 nb in S+U at 200 GeV and 0.54 nb for Pb+Pb at 158 GeV in the NA50 
spectrometer acceptance, 0.4 > xp > and —0.5 < cos 9 < 0.5 (to correct to the full 
angular range and 1 > xp > 0, multiply these cross sections by ~ 2.07). The experimental 
results in this range are 1.03 + 0.04 + 0.10 nb for S+U collisions § and 0.44 + 0.005 + 0.032 
nb for Pb+Pb reactions [[[]]. Interestingly, in the Pb system we find a Drell-Yan cross 
section (a/AB) DY = 35.2 pb while NA50 finds (a/AB) DY = 32.8 ± 0.9 ± 2.3 pb. Both 
the if) and Drell-Yan cross sections in Pb+Pb collisions are somewhat above the data, 
suggesting that the calculated rates at the NN level may be ~ 20 — 30% too large at 
158 GeV. Such an error is consistent with ambiguities in current pp data near that low 



energy [|TIJ. Note also that nuclear effects on the parton densities omitted in eqs. (4,5) 



can affect the total S and Pb cross sections at this level. 



4. COMPARISON TO OTHER WORK 

To compare more easily our results with those in Refs. |1]|§, we present the follow- 
ing technical points. Blaizot and Ollitrault [[| and Kharzeev in these proceedings have 
observed that one can describe the measured total cross sections in pA and S+U in 
the absence of comovers if one takes u^n = 6.2 mb. We also find this to be true - 
neglecting comovers and taking their value for the absorption cross section, we obtain 
(a/AB) t = 1.05 nb in S+U at 200 GeV and 0.62 nb for Pb+Pb at 158 GeV. On the other 
hand, our calculated _&r~dependent ratios do not agree with measurements, as shown in 
fig. 5 (see also ||). The importance of comovers in high Ep events follows from the cor- 
relation between Et, centrality and particle (comover) multiplicity. On the other hand, 
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Figure 5. Suppression due to nucleon absorption alone using a value ct^n = 7 mb as 
suggested by Kharzeev and Satz [2] does not describe S+U data (b). Both this value and 
our preferred 4.8 mb are consistent with pA data (a). 



comovers only contribute to the total S+U cross section at the ~ 18% level, because 
the impact-parameter integrated cross section is dominated by large b and the distinc- 
tion between central and peripheral interactions is more striking for the asymmetric S+U 
system. 

How well does our value ct^jv =4.8 mb agree with the pA data? We estimate the \ 2 P er 
degree of freedom to be 0.34 for = 0.62 mb and 1.3 for a^ N = 0.48 mb. Both values 
of x 2 /d.o.f. are close to unity and indicate reasonable fits. Given the paucity of the data 
and the crudeness of the models, we consider the agreement to data to be equivalent. 
Observe that including the measured S+U total cross section in extracting a^N from data 
as in [Q] introduces an inappropriate bias. 

Where does the L(Ejn) in NA50's plot come from and how does it approach the geomet- 
ric? To address these questions we estimate L(Et)- Following |T(J we neglect comovers 
and equate (3) to our result for the ijj—to- continuum ratio from (1) and (4-6). We then 
expand in powers of ct^n and equate the terms proportional to <j^n to find 

p L(E T ) = {2T AB y l J d 2 b P(E T , b) J d 2 s {[T A (s)} 2 T B (b - s) + T A (s) [T B (b - s)} 2 } (9) 

where 

T AB = J d 2 b P(E T , b) J d 2 s T A (s)T B (b - s) (10) 

and T A = J p A dz. This quantity is strictly meaningful for values of the mean free path 
Po&i/jN *C L, as are the results of [flOfl . Following jlUlOfl, we overlook the fact that this is 
not true for realistic parameter values, and obtain the result in fig. 6. We see that both 
our L and the average value of (|6|) obtained from P(ET,b) are in excellent accord with 
NA50 results in 0. To compare to NA50, we calculate L for p "~ * 0.132 fm~ 3 (a more 
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Figure 6. L(E T ) for S+U and Pb+Pb from (9,10) (solid curve) is in good agreement with 
the NA50 result (dashed curve) used to create fig. 1 (dashed curve). The dot-dashed 
curve shows the average impact parameter b obtained from P(EV,&)- The dashed line 
begins and ends at centroids of the lowest and highest bins measured. 



realistic value is 0.17 fm 3 ). More importantly, we see that saturation is nearly achieved 
in the symmetric Pb+Pb system even when the average h is rather large. 



5. i// SUPPRESSION 

To apply eqs. (4-6) to calculate the tp'-to-ip ratio as a function of E T , we must specify 
a^, N , oyjv, and oy co . Following Ref. |TT|, we use pp data to fix Ba^, N / Ba^ N = 0.02 
(this determines F^>). The value of u^n depends on whether the nascent ip' is a color 
singlet hadron or color octet cc as it traverses the nucleus. In the singlet case, one expects 
the absorption cross sections to scale with the square of the charmonium radius. Taking 
this ansatz and assuming that the ip' forms directly while radiative x decays account for 
40% of ip production, one expects cr^> ~ 2.1a,/, for interactions with either nucleons or 
comovers ||. For the octet case, we take a^jq ~ ct^n and fix oy co ~ 12 mb to fit the 
S+U data. In fig. 5a, we show that the singlet and octet extrapolations describe pA data 
equally well. 

Our predictions for Pb+Pb collisions are shown in fig. 6. In the octet model, the entire 
suppression of the ip'-to-ip ratio is due to comover interactions. In view of the schematic 
nature of our approximation to 5 co in eq. (8), we regard the agreement with data of singlet 
and octet extrapolations as equivalent. 
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Figure 7. Comover suppression of if>' compared to (a) NA38 and NA51 pA data [1, 9] and 
(b) NA38 S+U data [8] (filled points) and preliminary data [1]. 



6. SUMMARY 

In summary, the Pb data cannot be described by nucleon absorption alone. This 
is seen in the NA50 plot, fig. 1, and confirmed by our results. The saturation with L 
but not Ex suggests an additional density-dependent suppression mechanism. Earlier 
studies pointed out that additional suppression was already needed to describe the S+U 
results [BJ; recent data [|I]] support that conclusion (see, however, g). To study the S+U 
question further, experiments can search for saturation in super-central S+U reactions. 
More generally, while comover scattering explains the additional suppression, it is unlikely 
that that explanation is unique. SPS inverse-kinematics experiments (B < A) and AGS 



pA studies near the ip threshold |L5| can help pin down model uncertainties. 

In the context of our hadronic scenario, the success of the if) and if)' predictions imply 
comover densities perhaps exceeding ~ 1 fm -3 in Pb+Pb collisions. Similar densities are 
found in string and cascade models of heavy ion collisions. We must ask if such hadronic 
models are applicable at such densities — they must certainly be at the very edge of their 
applicability. Nevertheless, such models have proven very useful as predictive phenomeno- 
logical tools. The ability to predict hadronic backgrounds that we have demonstrated is 
a prerequisite for using charmonium to search for new phenomena at RHIC and LHC. 

We thank Claudie Gerschel and Michel Gonin for discussions of the NA38 and NA50 
data and Miklos Gyulassy, Geppetto McLerran and Ezra Pisarski for insightful comments. 
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